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Abstract 

Double parton fragmentation is a process in which a pair of partons produced in the short- 
distance process hadronize into the final state hadron. This process is important for quarkonium 
production when the transverse momentum is much greater than the quark mass. Resummation of 
logaritlims of the ratio of these two scales requires the evolution equations for double parton frag- 
mentation functions (DPFF). In this paper we complete the one- loop evaluation of the anomalous 
dimensions for the DPFF. We also consider possible mixing between the DPFF and single parton 
power suppressed gluon fragmentation and show that such effects are sub-leading. 
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The production of J/ip., T, and other quarkonium states in hadron colhders is an im- 
portant test of our understanding of perturbative QCD. Many early calculations used a 
color-singlet model (CSM) in which the heavy quark-antiquark pair {QQ) were assumed to 
be produced with quantum numbers identical to the final state hadron, e.g., ^5*^^^ for the 
J /ip and T. The CSM calculates the production cross section for producing a QQ pair with 
the right quantum numbers, multiplied by the quarkonium wavefunction. This model had 
to be rejected when measurements of prompt J /if) production at large transverse momen- 
tum at CDF revealed order of magnitude discrepancies between the CSM and experimental 
data. The modern theory of quarkonium production is based on Non-Relativistic QCD 
(NRQCD) In this theory, the QQ pair can be produced in a state with arbitrary quan- 
tum numbers. The nonperturbative transition of the heavy quark-antiquark pair produced 
in the short- distance process into the final state is governed by the QCD multipole expan- 
sion, and hence is an expansion in v, where v is the relative velocity of the heavy quark and 
antiquark. This theory allows for new color-octet production mechanisms that can natu- 
rally explain the size of the cross section at large (2| but fails to explain some features 
of the data, most notably the polarization of the J /if) at large For an up-to-date 

next-to-leading order analysis of quarkonium production with comparison to a wide variety 
of experimental results, see Refs. jj-l6[. 

At large p\_ the dominant mechanism for quarkonium production is fragmentation, in 
which we consider the probability to produce an "on-shell" parton which then has some 
probability to hadronize to the quarkonium of interest. To make sensible predictions for 
such processes one first needs to prove factorization, which, at least in this context of this 
paper, implies that at large the cross section can be written as the product of probabilities. 
The parton to hadron transition is encapsulated by a fragmentation function which obeys 
a DGLAP evolution equation that can be used to resum large logarithms of p^jviiQ that 
arise in the cross section. Within the framework of NRQCD, this resummation has been 
performed for single parton fragmentation contributions to quarkonium production 04^- 

There are many power suppressed corrections to this result, which naively one might think 
are numerically irrelevant. However, for quarkonium fragmentation there is a particular 
power correction, coming from double parton fragmentation, whose contribution can be large 
for moderate values of p^. While these corrections are down by rriQ/p^^ they are enhanced 
by a factor of compared to the single gluon fragmentation function. In fact, the double 
parton fragmentation operators are the unique power correction with this enhancement. 

For the case of quarkonium production, logarithms of p^/mq in double parton fragmen- 
tation have not been resummed. The double parton fragmentation function (DPFF) was 
first introduced in Refs. [l^, where it was evaluated in perturbation theory. It was 
shown to account for 80% of the next-to- lea ding order CSM cross section and furthermore 
yield longitudinally polarized J/ if). In Ref. [12] we derived factorization theorems for the 



DPFF using Soft Collinear Effective Theory (SCET) [13|, ll4|] and computed the evolution 
equations for the DPFF. In that paper, we focused on color-singlet and color-octet ^5*1 
production only. However, production in other spectroscopic channels such as ^S'q^'' and 
"^Pf^ are also important to collider production of quarkonia, at least at moderate trans- 
verse momentum 15, 1^. Furthermore, the DPFF can mix into power suppressed single 
parton fragmentation functions. The purpose of this paper is to complete the leading order 
calculation of the anomalous dimensions for the DPFF. 

The DPFF, D'^^{u,v,z), was defined in Ref. [12| in terms of the following SCET matrix 
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element: 



= 8 5{u[ — U2 + Lo'n — f — dudv 5{z ^ )6{v — 1 — z—^)6{u — z—^) 

J z cui — CO2 n'-p n'-p 

xzDf^^iu,v,z), 

where a = 1,8 indicates the color-structure of the operator, and F-*^^ G |{^', '^'ife, ^'7!/}, 
where 7^, = 'j'^ —n"^'^' /2 — n"^'i^' /2. The notation used here is the same as in Ref. [l2| and we 
refer the reader to that paper for details. The matching of DPFF operators onto NRQCD 
production operators is also discussed in Ref. At lowest order in v, the operators 

with Dirac structures Fi, F2, and F3 match onto ^^i, ^S'o, and ^Si NRQCD production 
operators, respectively. Fi corresponds to longitudinal polarization of the QQ pair, while 
F3 corresponds to transverse polarizations. P-wave operators can appear at higher order in 
the V expansion. In Ref. [12], only the Dirac structure F3 = |^'7X' was considered. It does 
not mix with the other Dirac structures under renormalization. We will see below that Fi 



and F2 mix among themselves. The variables z, u, and v are also defined in Ref. [12 



n'-p n'-p 



u'l — Uj'2 — W3 

z—^ = 1 + z-^ 
n'-p n'-p 

OJ'a ^ Oj'n 



u = z— — = 1 + 2;- 



n' -p n'-p 



Here, z corresponds to the hadron if's fraction of the QQ pair light-cone momentum, and 
u and V correspond to the fraction of the total QQ light-cone momentum carried by each of 
the heavy quarks in the QQ pair. Eq. ([T]) can be inverted to obtain 



D^^{u,v,z) = o / dijj'^dbj'2djjj''>diJi^b{{jj'Y — ijj'2 ^)^('^2 ^(^ — 1)) ^(1^4 — 



-u 



z 



The anomalous dimensions of the DPFFs are determined from the ultraviolet (UV) coun- 
terterms to Df^{u,v, z), which are calculated from the one-loop diagrams given in Fig. [H 
We work in Feynman gauge and use dimensional regularization to regulate UV divergences 
in loop diagrams. Here as in Ref. [l2|, SCET coUinear fields are represented by single lines 
and the double lines correspond to Wilson lines. There are additional diagrams that can 
be obtained by reflecting diagrams in Fig. [T] about the horizontal or vertical axis, or both 
axes. We refer to diagrams obtained by reflecting about the horizontal axis by adding a 
hat, e.g. C, diagrams obtained by reflecting about the vertical axis by adding a bar, e.g. C, 

and by doing both reflections by adding a hat and bar, e.g. C. Note that A = A, D = D, 
and E = E so these are not distinct diagrams. All other possible reflections lead to new 
diagrams. 

In Fig. [H we define the outgoing momentum of the quark and anti-quark on the left-hand 
side of the cut to be p'^ and P3, respectively, and we define the incoming momentum of 
the quark and anti-quark on the right-hand side of the cut to be p^ and P2, respectively. 
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FIG. 1: The one- loop diagrams required for computing the evolution of the DPFF. Not shown are 
the diagrams which are mirror images with respect to the horizontal and/or vertical axes. 



The large components of these momenta are expressed in terms of the momentum fractions 
defined in Ref. 12 : 



P 



X 

X 



n' ■ {pi +P2) 



(4) 



X- 



n ■ Pi 



X- 



n' ■ pi 



P 



Here P is the large light-cone momentum component of the final state QQ pair. 

Diagram A and its reflections have different color factors depending on whether or not 
the operator is color-singlet or color-octet. These color factors are where (5^^^ = Cp 

and = — Evaluating diagram A yields 



MX 



a, 



27r e V2P 



u 



u 



— zlx)b{i-v)0'^'' -e{\-u)^^Q[u-\) 



X 



X 



(5) 



On the right hand side the index a is not summed over. We define u = l — u, v = l— v, 
A = 1 — A, and ^ = 1 — ^, and the operator O*'" is 



(6) 



The diagram related to A by reflection about the vertical axis is obtained by making the 
replacements 



Mf = M'/{u^v,X^O 



(7) 
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The result of evaluating diagram B is 



Pis Cj. 



1 r zu 
- + ln — 



z/x)S{\-u)5{^-v)0''\ (8) 



Note that this diagram exhibits a rapidity divergence that has been regulated using the 
regulator of Refs. 171, ll8|. This is the origin of the factor of 1/t], where rj is the rapidity 
regulator, as well as the scale u. We will see below that the rapidity divergences cancel in 
the sum over all diagrams. The diagrams related to B by symmetry are obtained by making 
the following replacements: 



Mr 

B 
B 



M'j^l'iu ^ M, A ^ A), 



(9) 



Both color structures give the same result so the sum of all diagrams is 



^ia ^ia ^i,a ^i,a 
^ B B B 

as Cp 



4 , 
- + In 

7] 



4 4 
Z U 



UUVV 



pi 



+ 4 



(10) 



— ) 6{l-z/x)6{X-u)6{^-v)0''\ 



Diagram C evaluates to 



(11) 



a. 



a) 



27r 



1 r zu 
- + In — 

r] \M_r 



u9{u-\) 

6{X-u) - ^- — 

X[u — A)+ 



2F J X 



where again a is not summed over and the diagrams related to C by symmetry are obtained 
by making the substitutions in Eq. The result of summing diagram C and its reflections 
is 



a. 



C 



C 



c 



4 , f z^u 
- + In _ 

1] \uuvvP* 



271 e 

u e{x-u) u e{u-x) 

A (A — A (A — 



6iX~u)Si^-v) 



(12) 



<5(e - v[ 
6{X-u) 



2P 



5(l-^/x)0^'^ 



The virtual diagrams and their reflections do not lead to any mixing between singlet and 
octet operators or between Dirac structures 1 and 2. The l/rj poles and the corresponding 
logarithms of u cancel between diagrams B and C and their reflections when the operator is a 
color-singlet, but not when it is a color-octet. In the color-octet case the rapidity divergences 
cancel against rapidity divergences in the real emission graphs which we evaluate next. 

Real radiation comes from diagrams D — G and their reflections. The color factors for 
these diagrams are given by: 



ab 



(3u (3i8 

1^81 1^88 




ab 



Cf 

2Nr 
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The first index, a = 1 or 8, refers to tlie color state of the initial and final state quarks in 
the diagram and the second index, b, refers to the color-structure of the operator. The real 
radiation diagrams shown in Fig. [Devaluate to: 



Ml 



x'^ 1 — z/x 



X 



6{v--n6{u--\ 



X 



(13) 



271 \2PJ ec/y'^'^ 



x^ \ — zjx 



X 



X 



27r \2P) euv 



ab 



V 



+ In 



zu 

uP 

z 



5(1 - z/x) + 



ux 9{1 — z/x) 
\z (1 — z/x)j^ 



x5[v--n5[u--\]0 



M, 



G 



2n\2PJ euv 



A 



ab 



7] 



+ In 



X 

zv 
z 



(5(1 - zjx) + 



ux Q{\ — zjx) 
\z (1 — zjx^j^ 



X 



^b[v- -i\b[u- -\\0 



Z V 



where h is not summed over. The reflections oi D — G are given by the substitutions in 
Eq. ([H]). The real diagrams mix color-singlet and color-octet operators, and diagrams D and 
E also mix the Dirac structures. 

The DPFF is renormalized multiplicatively, 

^QQ(bare)^^^ t^, ^) = ^ dv! dv'^^Z.a.biu, u', v, v', z/z', li)Df{u\ v', z', /i) , (14) 

and thus obeys the renormalization group equation 



^^^Df^?{u"y',z",^,) 



du'dv'—^,a,bW\ u\ v'\ v\ z"/z', fi)Df;^iu\ v\ z\ /i) , (15) 



where the anomalous dimension is given by 

'yia,jb{u,u',v,v',z/z',fi) = 



(16) 



dz" d 
du"dv"—Zr^j^^{u, u", v, v", z/z", n)n—Zkc,jb{u", u', v", v', z" / z', n) . 

The indices i and j label the Dirac structure and the indices a and b refer to color (a, b = 
lor 8). The tree-level matrix element of the DPFF using partonic states with momenta 
labelled by {x, A,^) is given by 



D'',^{u,v,z) 



5(1 - z/x)6{X - u)6{^ - v)0^''' 



' j,b V"' - ' -/ \2P J 

Finally, we need to include the wave function renormalization, 

asCp 



1 



47re ' 



(17) 



(18) 
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after which we find that the anomalous dimensions are: 

asCp 



711,11 



7l8,21 



^5(1 - z/z') (sSiu - u')S{v - v') 



(19) 



u 



+6{v-v')\e{u' -u)- 



+5{u-u'){e{v'-v)- 



Xu' -u)+ 

1 

{v' - v)+ 



+ 1 
+ 1 



II 

+ e{u-u')- 

u' 

+ e{v-v')-, 

v' 



[u — u') 

1 

{v - v')^ 



- + 1 

h 

+ 1 



721,21 = 7ll,ll ) 

718,11 = -^e{i - z/z') 



(20) 



uv'-\-vu' ^1 — z/z' 



(1 - z/z')+ 2zlz' 



1 _ z _ 



U'V 



+ 



uv' + vu' 1 — zl z' 

+ ' 



{l-zlz')+ 2z/z' 
uv' + vu' \ — zj z' 

(1 - zjz')^ ^ 2z/z' 
uv' + vu' \ — zjz' 



XI- z/z') + 2z/z' 



-^,6{v - -v')5{u - -u') 
u'v' z' z' 

\ z z 

5{v v')5{u u') 

u'v' ^ z' ' ^ z' ' 

-^,S{v - -v')5{u - -u') 
u'v' z' z' 



(21) 



uv' + vu' 1 — z/ z' 



, , , , — 5iv v')5iu 

_{l-z/z')+ 2z/z' J u'v' ^ z' ' ^ 



+ 



uv' + vu! 1 — z/ z' 



+ 



{l-z/z')+ 2z/z' 
uv' + vu' 1 — z/ z' 

(1 - z/z')+ ~ 2z/z' 
ur' + ru' 1 - z/z' 

Xl-z/z')+ ~ 2z/z' 



■^^5{v - -v')5{u - ^^u') 



u'v 



—5{v--v)6{u-- 



u 



u'v 



—6{v - -v')6{u - -u) \ , 



u'v 



728,11 — 7l8,21, 
728,21 = 7l8,ll- 

In addition, 



7a,i8 



2N, 



■7i8,il 



(22) 

(23) 
(24) 



(25) 
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Furthermore, 



718,18 = -?,'^5{u-u')5{v-v')5{l-z/z') 

TT 



(26) 



u 



6{u' -u)—\ — ^ h 1 I + (u ^ u,u ^ u'^ 



u' \ {u' — 



+ -Jtt6{u-u')5{1-^) 



TT 2N, 
TT \Z' 



z' 



[V — V) — 



V \[V — V) 



+ 1 ) + (f ^ 17, f ' ^ f ') 



uv' + vu' \ — zl z' 



-^b{v- -v)b{u - -u) 
u'v' z' z' 



+ 



2N^ [{l-z/z')+ ' 2z/z' 

-2 r UW' + Wm' 1 - z/z'] 1 2 ,^ 2 n 

7 H r — v')6(u u') 

1 - z z')+ 2z z' J u'v' ^ z' ' ^ z' ' 



2Nr 



1 
1 



uv' + vu' \ — zl z' 
+ ' 



(1 - zjz')^ 2z/z' 

uv' + vu' \ — zj z' 
Xl-z/z')+ ^ 2z/z' 



u'v' ^ z'^' ^ z' ' 
^,5{v-^-v')5{u-^u'))e{l-z/z% 



7l8,28 — I — 7 

77 V^;' 



728,18 
728,28 



a, / 2 \ 2 { Ni -2 



2Nr 



uv +vu 1-z/z] 1 ^, z z 

, TT^H T-r- —r-Av -v)b{u -u) (27) 

1 - z z')^ 2z z' J u'v' ^ z' ' ^ z' ' ^ ' 



+ 



+ 



iV2-2 



2Nr, 



uv' + vu' 1 — zl z' 



(1 - z/z')^ 2z/z' 



uv' + vu' \ — zj z' 



1 



{l-z/z')+ 2z/z' 

uv' + vu' \ — zj z' 
{l-z/z')+ ~ 2z/z' 



z 

U'V ■ Z' 

z 



—5{v--v')5{u--u') 



-u 



uv 



U'V 



Si^^-^')Si^u--u')\e{l-zlz') 



7l8,28 
7l8,18- 



(28) 
(29) 



In addition to the above mixing matrix we must also consider the mixing with power 
suppressed single parton fragmentation functions. As emphasized in the introduction, while 
the direct contribution of such operators will not receive any 1/v enhancement, in principle 
they can still affect the runnings of the DPFFs. For instance we may consider the mixing into 
the power suppressed single gluon fragmentation function. Such mixing can in general arise 
from diagrams such as those in Fig. |2^. However, to mix from the single gluon to the double 
quark fragmentation function would require diagrams such as those in Fig. |2]d, which is sub- 
leading in a{p_[_). Thus at lowest order the anomalous dimension matrix will contain only 
one off-diagonal element that will not generate any logs proportional to the DPFF. Then, 
given that the single gluon power suppressed fragmentation function is not enhanced we may 
ignore any mixing between double and power suppressed single fragmentation functions. 

To summarize, we have extended the calculation of the anomalous dimensions of the 



DPFF in Ref. [12| to include all possible Dirac structures. This gives the complete set of 



anomalous dimensions needed at this order in p±/mQ. The renormalization group equations 
presented here and in Ref. [12|] can be evolved from the scale p± to rriQ to resum logarithms 
oip±/'mQ in double parton fragmentation contributions to quarkonium production. When 
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r 



r 



'U"0"0"0"(J 



o 




a 



b 



FIG. 2: Possible mixings between single and double parton fragmentation functions. 



this is completed, it will be interesting to see if the resummed cross sections can improve 
our understanding of quarkonium production. 
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